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^ ■ Abstract 



Magnetic resonance in ordered state is shown to be the direct method for distinguishing the 
orbital ground state. The example of perovskite titanates, particularly, LaTiOs and YTiOs, is 



D ' considered. External magnetic field resonance spectra of these crystals reveal glaring qualitative 



dependence on assumed orbital state of the compounds: orbital liquid or static orbital structure. 



j^ , Theoretical basis for using the method as an orbital state probe is grounded. 
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FIG. 1: RTiOs crystal structure. Pnma (a, b, c) and pseudocubic (x, y, z) axes are shown. 
Numbers denote Ti sublattices. 

Recently, a lot of spin and orbital phases and their phase transitions attract intent inves- 
tigation due to interplay of these degrees of freedom, particularly in transition-metal (TM) 
oxides, pj Fundamental physical properties revealed in such systems are still a subject for 
a discussion. Among the phenomena which attract the most attention there is a superex- 
change interactions driven rich spin-orbital quantum phase diagram proposed to exist in 
perovskite titanates and vanadates. |^, y, |j] 

Whether orbital liquid state present in real compound or not — is the question which 
have given rise to hot debates especially concerning the simpler system — RTiO^, R is 
rare-earth element or Y. 

In wide temperature range for different R, titanates are known to possess orthorhombic 
crystal structure |5|, |6| (which is often called "quasi-cubic") with Pnma space group (see 
Fig. [1]). GdFeOs-type distortions (Ti^-distortions), which are present in these crystals, are 
believed to control magnetic structure and properties of the compounds through the influence 



on their orbital ground state. 



U,mM 



The discussion of titanates orbital ground state was started by Khaliullin and Maekawa 
who proposed the superexchange model with dynamical quenching of local orbital moments 
in a simple cubic lattice of Ti^+ ions with no static orbital order but with fixed magnetic 
arrangement — in contrast to usual Goodenough-Kanamori picture. This approach based 
on the Kugel-Khomskii model [9|] perfectly explains the unusual reduction of Ti'^"'" spin 10|, 



Ill ] and anomalous isotropic spin-wave spectrum [12] found experimentally, 

7, 



3ut contradicts 
I, y, bj and 



to NMR [IJ and XAS UM experiments as well as some crystal-field 

nn 

density-functional |17l . Ilq calculations. However, it is supported by recent Raman scattering 



experiments [19| revealing orbital excitations and thus making some orbital fluctuations 
possible. 

In this context La and Y titanates are believed to be of special interest for investigators 
as these two ions stand at the opposite ends of rare-earths and Y series with different ionic 
radii. Ti^-distortions in YTiOa with "smaller" Y are much greater of these in LaTiOa, 
which has the biggest radius of R-ion in the whole series. One more feature which makes 
these two crystals more interesting is their magnetic ground states: lanthanum titanate is 
antiferromagnetic |6| with strong isotropic superexchange J of about 15.5 meV, [l^ whereas 
yttrium titanate is ferromagnetic with J ^ —2.75 meV. ll|] Moreover, almost isotropic 
superexchange couplings both in LaTiOa and its sister compound are in contrast with the 
situation in manganites and cuprates which have pronounced superexchange anisotropy. 

Keeping in mind a huge collection of all above mentioned facts, one still can not find 
a compromise between them and can not answer what ground state is actually realized in 
titanates. That indicates a puzzling situation: on the one hand the system is already studied 
in many aspects, on the other hand some of its properties are not clear yet. To meet this 
challenge, the definite model of a titan perovskite oxide is developed in the current paper 
on the particular example of Y and La compounds 

We start from the crystal-field Hamiltonian with explicit electron-lattice interaction, that 
is vibronic Hamiltonian (VH): 
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Hlin + Hqq + Hr = 

[Ve {QeXe + QeX,) + Vt {Q^X^ + Q,X, + Q^X^)] + 
Va (QI + QI + Qt) Xai + V,{{2Ql -Ql- Ql) Xe + V3 {qI - Ql) X,} + 

Vc [QyQzX^ + QxQzXn + QxQyX(^)] 
Ve"" {QfXg + gf X,) + 

Vt"" { {Qt + QI2) Xc + (Q^i + Q^,,) X, + {Ql, + Ql,) X^} 

Here Qr and Qr (r = 9, e, ^, rj, (, x, y, z) are symmetrized shifts of oxygen and i?-ions which 
are nearest and next-nearest Ti'^"^ neighbours correspondingly. These shifts are obtained 



(1) 



y and YTiOs. m 



from accurate crystal structure data for LaTiOa [6j and YTiOs. [5| Xr are symmetric orbital 
operators, acting on the 3d-t2g triplet, and Va {a = e,t,a,b,c) and V^ (a = e,t) are 
electron-lattice coupling constants. |20|] The first two square brackets, which are Hun and 
Hqq, represent Ti'^+ 3d-t2g electron interactions with nearest oxygen linear and quadratic 
symmetric shifts correspondingly, whereas Hji is responsible for this electron and -R-ions 
shifts coupling. 



This Hamiltonian with ab initio Va (calculated with GAMESS package 2l|, |2^) gives 
Tound state with static orbital structure of the form predicted by Mochizuki and Imada, [7|, 
I that is the following orbital functions {3d-t2g cubic basis set): '?/'i,2(La) = ip3^4^(La) ~ 
l/V3{—^ — ri±(), which is almost trigonal, in LaTiOs and '?/'i,2(Y) ~ 1/v^ (=F^ + C); 
V^3,4(Y) ^ l/-\/2 (=i=?7 + (^) — in YTiOs. Here sign alternation reads sites 1(2) - the up- 
per sign, and 3(4) - the lower one. 

We have irrefutable argument for reproducing results of previous LDA-based [ij, [l7|, [l8| 
or point-charges [Tl, M, UM Il6 1 investigations by using ([T]). All these studies either were 

71, Isj (giving unhkely structure parameters) or did not reveal 

n n n ^^ 



based on oversimplified model 

the mechanisms of the ground state formation, p^. 
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18| The present calculation 



based on ([T]) showed that Hji plays crucial role in the formation of solitary orbital singlet 
with its segregation at ~ 0.20 eV in LaTiOs and ~ 0.15 eV in YTiOs. A half of these gaps 
is produced by the R-ion crystal field. The influence of the remainder of the crystal on Ti^+ 
orbital state is negligible. 



we then exploit common Kugel- 
9| Thus, arriving to isotropic superex- 



Using the low-energy spectrum obtained from 
Khomskii (KK) method within the Hubbard model, [l 
change, one then can perform Moriya's approach [23i] for treating antisymmetric terms ^] 
of the effective S-^ spin-Hamiltonian (ESH) introduced below: 



H,ff = J,, (S, ■ S,) + D,, (S, X SO + S^A^fS^ + H^gfS', 



^i 1 



(2) 



where Jij stands for isotropic superexchange between z-th and j-th magnetic ions, Djj is 
Dzyaloshinskiy-Moriya vector, Aj^ is symmetric anisotropy tensors, gj is g-factor and H 
represents external magnetic field. 

The Hubbard model parameters: energy of electron hopping from the m-th orbital on 
one site to the n-th on the neighboring one in strictly cubic system tmn, on-site Coulumb 
repulsion of a pair of electrons U and intra-atomic electronic exchange interaction Jh — 




FIG. 2: Schematic magnetic structures of LaTiOa (up) and YTiOs (bottom). Gray arrows denote 
magnetization projections (Pnma). 



TABLE I: Single crystal g- factors of LaTiOs and YTiOs. The Pnma coordinate system is used. 



R 



La 



Y 



' L93 0.03 -0.05 



0.03 L89 -0.04 
-0.05 -0.04 L93 



V 



L93 -0.04 -0.04 
-0.04 L87 0.02 
-0.04 0.02 L96 



are taken from LDA-based calculation. [17[ For LaTiOs tmm ~ 0.25 eV, t^^n ~ 0.12 eV, 
U ^ 3.20 eV, Jh ~ 0.61 eV. For YTiOa t^m ~ 0.27 eV, t„^„ ^ 0.08 eV, U ^ 3.45 eV, 
Jh ~ 0.62 eV. Indexes in t^n denote 3d-t2g orbitals, namely ^, rj, (. All interactions in ESH, 
which follow from these parameters, are listed in Tables [H [TTl 

This calculation of superexchange couplings, although it is not something new, is to be 
reproduced because of extreme sensitivity of these parameters to the orbital state. This well- 



known feature o 
studies [Tl, l8|. 
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the KK treatment should be considered as (in comparison with previous 
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25j) we have obtained new ipilLa) and ipiO^) 



1...4. 



By using the Hamiltonian ([2]) one can obtain magnetic ground state as well as magnetic 
excitations in both compounds. The magnetic structure type for both crystals is {A^, Fy, Gz) 
with major G-component in LaTiOa and F — in YTiOa (see Fig. [2]) — in excellent agreement 
with neutron scattering experiments. 10|, lllj Here we should emphasize two special features 
in Jij which are crucial for obtaining correct magnetic structure. First is considering Hund's 
coupling, which is Jh/U smaller then "common" superexchange, [26!] proportional to t^n/U 
for titanates. The second is explicit introduction of the Ti-O-Ti bond angle (</)) -dependence 



TABLE II: Magnetic interactions (nieV) in LaTiOa and YTiOa. The Pnma notations are used. 



R Ji2 Jl2 



D 



12 



D 



13 



^12 



U3 



La 13.21 16.12 



Y -2.77 -2.72 



V 



0.58 



-2.80 

0.60 


-0.21 



V 




-0.36 
-0.17 


-0.38 
-0.11 



-0.07 0.06 

-0.08 

0.06 -0.34 

-0.21 0.06 

-0.03 

0.06 -0.06 



V 



-0.02 

-0.03 -0.02 

-0.02 -0.01 

-0.69 

-0.10 -0.09 

-0.09 -0.16 



of tmn- 271] These two give for a pair of Ti'^"*" ions along the b axis {Pnma): 

Ji2 = a\b 



Jh 
U 



-c{d+ -^) sinV, 



(3) 



where a = 3.54(2.61), b = 6.37(1.66), c = 43.30(14.70) and d = 0.84(0.92) are numeric 
coefficients for La(Y) compound. These coefficients depend on the particular Ti'^+ orbital 
state, on t^n and U. Ti-O-Ti bond angle if for lanthanum titanate is about 153° and 142° for 
YTiOs. Without Jh one can not obtain experimentally observed magnetic structure in both 
compounds simultaneously as Schmitz et al. [l5|, [l6| couldn't. And in strictly cubic system 
(without sin^(y9-dependence of Jij) one is not able to obtain the correct orbital structure as 
KhaliuUin and Maekawa couldn't [2] for LaTiOs. 

It is important to mention that both anisotropic terms of the ESH, namely 
Dzyaloshinskiy-Moriya interaction and symmetric anisotropy, should be considered, oth- 
e™=eihe=tat.c .nagnet.c orie. couldn't e.i=t. fl Th. i. not always kept by inve^ga- 

tors. \m\m 



IDC tlJ 
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Actually, the above result is not unique as its different parts were obtained by several 
[3, Isl, ll5|, [la, llTI, [l8|, |25|] It is reproduced to illustrate the realistic model with 



authors 



reasonable parameters. But this model reveals particular mechanisms of titanates orbital 
and magnetic ground state formation that was not performed before. 

Now we turn to the core idea of the paper. That is drastic dependence of the particular 
kind of magnetic excitations on orbital ground state. 

We consider magnetic excitations, namely spin waves (SW) and antiferromag- 
netic/ferromagnetic resonance (AFMR/FMR) spectra, in LaTiOa and YTiOs exploiting 
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FIG. 3: (Color online) Spin wave dispersions for LaTiOs (up) and YTiOa (bottom). Circles denote 



experimental data, taken from Ref. 



12 for LaTiOs and from Ref. 11 for YTiO^. Simulated curves 



are black for static orbital order and dashed red for orbital liquid. 



linear approximation. Simulations (see Figs. [3], HI) turned out to be surprising when com- 
paring the real static orbital structure with hypothetic orbital liquid state. The latter was 
simulated by averaging the superexchange parameters through low-energy part of the spec- 
trum ([1]). 

There is an unexpected feature: spin-wave spectra exhibit almost no changes in two 
different orbital states for both crystals. Little discrepancy in those simulations can be 
easily removed by the slight fitting of the Hubbard model parameters, which fitting is re ally 



possible within these parameters calculation discrepancies in different approaches. [17l.ll8U29| 
At the same time there is drastic change in AFMR/FMR spectra for La and Y compounds. 
Orbital liquid shows up here in two ways: rising anisotropy in LaTiOa and suppressing it in 
YTiOs. This produces handbook curves of AFMR and FMR field spectra, but with different 
anisotropic behavior. 

This result is due to killing electronic distribution anisotropy factor in forming the mag- 
netic interactions by the orbital liquid. Thus, the lattice becomes the only source of magnetic 
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FIG. 4: (Color online) Field dependencies of AFMR frequencies in LaTiOs (up) and YTiOs (bot- 
tom). H is parallel to c for LaTiOs and b for YTiOa (Pnma). These directions are close to easy 
axes. Curves are black for static orbital order and dashed red for orbital liquid. 

anisotropy. That is why the latter appears to be different in the compounds under consid- 
eration — quite expectable for La and Y crystals. 

Unfortunately, no any attention was paid to such powerful and sensitive method of mag- 
netic structure and magnetic couplings investigation as antiferromagnetic/ferromagnetic res- 
onance so far. The only attempt to observe electron spin resonance (ESR) below the mag- 
netic transition temperatures (this is AFMR or FMR) and above them (EPR) was made by 
S. Okubo et al. [3Q| But, firstly, the AFMR signals for LaTiOs were not obtained at all, that 
might be because of frequency limits of the equipment used. And, secondly, this investiga- 
tion was made only for powder samples and thus all direction-dependent effects were wiped 
out. We can't but mention that if powder samples are used one neither will observe such an 
interesting field spectra as represented in Fig. HI nor the strong orbital state dependence of 



8 



these spectra can be obtained. 

Finally, we argue that up-to-day, magnetic resonance is the ultimate method for distin- 
guishing between static orbital order and orbital liquid. This method is a referee between 
opposite orbital states. Particularly, AFMR/FMR experiments with single crystals should 
put a dot at the end of the discussion what is the real orbital ground state in titanates — the 
remarkable system with strong entanglement of lattice, spin and orbital degrees of freedom. 

We thank I.V. Solovyev for fruitful discussions. This work was partially supported by the 
"Dynasty" foundation, by the CRDF REC-005, by the Russian Federal Scientific Program 
"The development of scientific potential of the high school" and by the RFBR grant no. 08- 
02-00200. 
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